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ON  CONTINUOUS  DISTRIBUTIONS  OF  DISLOCATIONS  IN  NONLOCAL  ELASTICITY 


A.  Cemal  Eringen 
Princeton  University 
Princeton,  NJ  08544 


ABSTRACT 

A  linear,  nonlocal  continuum  theory  of  dislocations  is  developed. 
The  field  equations  are  given  for  the  dislocation  density  and  the  stress 
fields  due  to  continuous  distribution  of  dislocations.  Green's  functions 
are  obtained  for  two  and  three-dimensional  media  and  an  integral  formula 
is  given  for  line  distribution  of  dislocations  generalizing  Peach-Koehler 
formula  of  the  classical  (local)  theory.  Unlike  the  classical  theory,  no 
stress  singularities  occur  so  that  self-stress  and  energies  of  disloca¬ 
tion  loops  can  be  calculated  involving  no  divergences.  Exact  solutions 
given  for  the  line  and  circular  distributions  of  dislocations  verify  these 
expectations. 
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1.  INTRODUCTION 


In  classical  (local)  theory,  the  displacement  and  stress  fields 
due  to  a  continuous  distribution  of  dislocations  can  be  calculated  by 
means  of  various  surface  and  volume  integrals  once  the  Green's  functions 
for  the  dislocation  density  is  known.  To  obtain  the  Greens  function  one 
solves  partial  differential  equations  f^r  the  stress  functions.  For 
example,  the  stress  due  to  a  line  distribution  of  dislocation  is  given 
by  the  Peach-Koehler  formula.  One  of  the  basic  difficulties  of  this 
theory  is  tha*  the  self-stress  and  energies  of  dislocation  loops  possess 
mathematical  singularities  so  that  calculations  will  have  to  be  cut-off 
near  the  lines  of  dislocation  or  the  core  region. 

In  several  previous  papers  (cf.  1-4),  I  have  shown  that  the 
solutions  for  the  single  dislocation  involve  no  stress  or  energy  singu¬ 
larities  at  the  core  regions.  Moreover,  calculated  theoretical  strengths 
of  solids  and  dispersion  curves  for  plane  waves  agree  quite  well  with 
those  known  from  the  atomic  lattice  dynamics  and/or  experiments.  There¬ 
fore,  it  is  expected  that  the  nonlocal  theory  will  eliminate  the  classical 
singularities  for  the  self  stress  and  energies  of  dislocation  loops.  The 

A 

raison  d'etre  of  the  present  paper  stems  from  the  need  to  develop  a 
theory  of  continuous  distribution  of  dislocations  based  on  the  nonlocal 
elasticity,  that  can  hopefully  predict  the  physical  phenomena  in  the 
microscopic  and  atomic  scales  where  classical  theory  fails  to  apply. 

In  Section  2,  I  summarize  the  basic  equations  of  the  nonlocal 
theory  of  linear  isotropic  elastic  solids.  In  Section  3,  I  develop 


the  field  equations  for  the  continuous  distribution  of  dislocations. 

Green  tensors  for  the  Beltrami  and  Airy  stress  functions  for  solids  of 
infinite  extends  are  obtained  in  Section  4.  The  stress  fields  are  given 
by  a  volume  integral  generalizing  celebrated  Peach-Koehler  formula  of  the 
classical  theory.  In  Section  5,  I  derive  explicit  expressions  of  the 
stress  fields  for  the  uniform  distributions  of  screw  dislocations  along 
a  straight  line  segment  and  along  a  circular  loop.  The  exact  formulas 
obtained  for  these  cases  contain  no  singularities,  justifying  our  expec¬ 
tation.  Section  6  contains  calculations  of  stress  fields  and  a  discussion 
of  the  reduction  of  yield  stress  with  the  dislocation  pile-up. 

The  simplicity  and  the  aesthetics  of  these  results,  I  believe, 
justifiably  indicate  the  power  and  the  potential  of  the  nonlocal  theory  in 
the  treatment  of  the  physical  phenomena  with  characteristic  lengths  in  the 
microscopic  and  atomic  scales. 

2.  BASIC  EQUATIONS 

The  linear  theory  of  nonlocal  elasticity  is  based  on  Cauchy's 
equations  of  motion 

Vk  *  °(h  -  V  *  0 

and  the  integral  constitutive  equations 

(2.2)  tkA(x.t)  -  ‘k£nn (*’-*)  em(x')  dv(x') 

V 

where  t^,  p  ,  f^  ,  u^  and  e are  respectively,  the  stress 
tensor,  the  mass  density,  the  body  force  density,  the  displacement  vector 
and  the  linear  strain  tensor  defined  by 


In  (2.2),  c,  „  is  a  function  of  the  vector  x'-x  and  the 
integral  is  over  the  volume  of  the  entire  body.  Consequently,  the  stress 
at  a  point  x  depends  on  the  strain  at  all  points  x'  of  the  body. 

Throughout  this  paper,  we  employ  rectangular  coordinates  x^  , 
k=l,2,3,  and  use  the  usual  summation  convention  on  repeated  indices.  Also 
a  superposed  dot  indicates  the  time  rate  and  an  index  followed  by  a  comma 
partial  differentiation  with  respect  to  xk  ,  e.g. 

auk  duk 

\  s  W  ’  uk,*  =  ^ 


The  kernel  possesses  certain  symmetry  regulations  and  it  depends 

on  a  length  scale.  For  isotropic  solids  (2.2)  takes  the  simple  form 


(2.4) 


jQ(|x'-x|)  Ok£(x')  dv(x') 
V 


where 


o 

kg. 


is  the  classical  (local)  stress  tensor  given  by  the  Hooke's  law 


(2.5) 


k£ 


=  X  ®rr6k£  4  2lJ  ek£ 


and  a  is  a  function  of  the  distance  |x'-x|  .  It  also  depends  on  a 
length  scale  c  that  may  be  taken  to  be  proportional  to  an  internal  charac¬ 
teristic  length  a 


(2.6) 


e 


e0 


a 


where  eQ  1$  a  non-dimensional  material  property  which  may  be  determined 
by  one  experiment  or  comparison  with  calculations  based  on  lattice  dynamics3 


The  internal  characteristic  length  a  may  be  taken  as  the  lattice  parameter 
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for  single  crystals,  granular  distance  for  amorphous  materials,  and  the 
average  distance  for  fiber  composites.  As  t-+  0  ,  t^  -*•  and  (2.4) 

reduces  to  Hooke's  law  t^  =  0^  .  Thus,  a(|x'-x|)  is  a  Dirac  delta 
sequence. 

3  5 

In  several  previous  papers  ’  ,  I  have  discussed  the  properties  of 
a(]x'-x])  and  gave  representations  which  lead  to  excellent  agreement 
with  known  atomic  calculations  on  dispersions  of  waves in  the  entire 


Brillouin  zone  and  on  theoretical  strengths  of  solids, 
the  two-dimensional  case,  an  appropriate  kernel  is 


For  example,  for 


(2.7) 


oi(|x|,e)  =  (2Tre2)"1  ^(r-x-x/e) 


which  satisfies  the  equation 


(2.8) 


(l-eV)a  =  6  ( |  x '  -  x  | ) 


vanishing  at  infinity.  In  fact,  for  the  infinite  solid,  it  can  be  shown 
that  a  is  the  Green's  function  satisfying  (2.8)  in  three-dimensions  also. 
Using  (2.8)  in  (2.4),  we  obtain 

(2.9)  Cl-cV}^  ■  V 

By  means  of  (2.1)  and  (2.9),  we  then  find  that 

(2.10)  W*U)\,U2  *  V  *Mk  *  (1  -  e-V)(pft  -  p»t)  .  0 


These  are  the  partial  differential  equations  for  ,  replacing  Navier's 
equations  of  classical  elasticity.  For  the  static  case  and  vanishing  body 
forces,  this  reduces  classical  Navier's  equation 


S 


(2.11)  (XnO  ♦  u  ui>kk  =  o 


However,  note  that  the  stress  field  is  determined  by  solving  (2.9) 
under  appropriate  boundary  conditions. 


3.  CONTINUOUS  DISTRIBUTION  OF  DISLOCATIONS 

Continuous  distribution  of  dislocations  is  envisaged  as  follows: 

A  small  neighborhood  n(x)  of  x  in  a  distorted  body  of  volume  ,  may 

be  relaxed  to  a  small  neighborhood  N'(X)  of  the  image  of  X  of  x  ,  in 

an  undistorted  (or  a  natural)  configuration  V  ,  by  releasing  constraints 
exerted  to  n(x)  by  the  rest  of  the  body.  A  line  element  dx  at  xen(x) 
can  be  expressed  in  terms  of  its  image  dXCN(X)  by 


(3.1) 


dx  =  A  dx 


where  A(X)  is  called  the  elastic  distortion.  It  is  assumed  that  A(X) 
is  continuously  differentiable  and  possesses  uniqie  inverse  so  that 


(3.2) 


-1 

dX  =  A  dx 


f 
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Consider  a  smooth  surface  S  in  V  bounded  by  a  closed  curve 
true  Burger's  vector  b  of  the  dislocations  piercing  through 
fined  by 


(3.3) 


» 

-i 

(|  dX  =  ( 

)  A  dx  = 

c  c  s 

a  n  ds 


where  n  is  the  unit  normal  to  5  ,  the  positive  sense  of  C 
counter-clockwise  when  sighting  along  n  .  Here  a  is  called 
dislocation  density 


-1 

(3.4)  a  =  curl  A  or 


For  small  distortions,  we  can  write 


(5-5)  \t  '  5U  *  “k£  ' 


so  that 


a 


jk 


Ckmn 


-1 

A. 

jn.m 


-1 

=  6ki  '  ak£ 


(3.0)  a..  =  e.  a. 

jk  kmn  jm,n 

From  this,  it  follows  that 

(3'7)  ajk.k  “  0 


C  •  The 
S  is  de¬ 


being 
the  true 


The  linear  strain  tensor  e,  0  and  rotation  tensor  r,  p  are  given  by 


(3.8) 


ek£  =  2  (ak£  +  a£k5  ’ 


rk£  =  2  (ak£  "  a£,k^ 


The  strain  incompatibility  is  expressed  by 


(3.9) 


eijk  E£im  ein,jm  "  ^ki 


where  is  called  the  incompatibility  tensor  and 


is  given  by 


(3.10) 


r'k£  ~  2  ^ekmnan£,m  +  E£ninank , n? 


All  these  results  are  well-known  in  classical  theory  (cf.  Ref.  [7]). 

In  nonlocal  elasticity,  the  strain  tensor  can  be  solved  by  using 
(2.9)  and  C.S) 


(3.  11) 


ek£  =  2y  (1  ’  e  V  )(tk£  _  1+v  trr<Sk£-) 


where  v  =  A/2(A+p)  is  the  Poisson's  ratio.  Substituting  (3.11)  into 
(3.9),  we  obtain 


(3.12)  a  -eV)lv\t  .  ^  (trr_u  -  V2trr«ke)l  -  2u  nkt 


These  equations  must  be  solved  under  the  conditions  of  equilibrium 


(3.13) 


0 


8 


Following  Kroner's  approach  ,  modifying  the  Beltrami  solution  of  (3.13), 
we  take 


t3,14)  tk£/2u  =  7  Xk£  *  1-v  *-Xrr,k£  "  7  xrr  6k^} 


where  the  symmetric  stress  function  X^£  is  subject  to 


(3.15) 


Xk£,£  =  ° 


Substituting  (3.14)  into  (3.12),  we  obtain 


(3.16) 


(1  -  sV)  v\u  ■  v 


Thus,  given  the  dislocation  density  function  a^  ,  through  (3.10),  we 
calculate  .  The  solution  of  (3.16)  gives  and  (3.14)  the 

stress  field. 

Equation  (3.16)  is  singularly  perturbed  and  as  expected  in  the 
limit  e-*"0,  (3.16)  reduces  to  the  classical  equation  for  • 

To  obtain  the  solution  of  (3.16),  we  must  find  the  Green's  Tensor 

Gk£mn(*,P  which  satisfies 


(3.17) 


(1  -  e‘V“)  7  G, 


6 (x-C)  <5,  „<$ 

~  k£  ran 


The  solution  of  (3.16)  is  then  given  by 


(3.18) 


Gk£mn(?^ 


subject  to  supplementary  conditions  (3.15). 


4.  GREEN'S  TENSORS  AND  STRESS  FIELDS 


Here  we  determine  Green's  tensors  for  two  and  three-dimensional 
bodies  of  infinite  extends. 

(i)  Three-Dimensional  Infinite  Space 

2 

The  operator  V  is  invariant  under  rotations  of  coordinates.  For 
the  infinite  space,  we  look  for  a  solution  of  (3.17)  which  depends  on 
| x-C [  only,  i.e.. 


(4.1)  (1  -  e2V2)  V4G  =  6 (x  -  £) 

2  2  4 

Since  the  operators  1-e  V  and  V  are  commutative,  we  set 


(4.2) 


(1  -  eV)  G  =  H  , 

V4H  =  6(jx-C|) 


For  the  infinite  space,  H  is  given  by 


(4.3) 


H  =  -  j  x  -  Cl /8tt 


In  spherical  coordinates  using 


1  ,  2  d  . 

T  (r  d75 

r 


we  obtain  for  G  : 


' 
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e2G  |x-£| 

(4.4)  G(|x-||)  =  47]p|7  exp(-|x-C|/c)  -  ,  *  *  0 

(4.5)  G(|x-||)  =  -  .  e  =  0 


Here  is  an  arbitrary  constant  which  may  be  chosen  GQ  =  1  to  render 


t^  regular  at  x *  £  .  The  solution  of  (3.16)  for  the  infinite  solid 
is  given  by 


(4.6) 


xk£ 


r 

G('Hl)T1kt(S)  dv(^ 


v 


which  satisfies  the  conditions  (3.15)  on  account  of  (3.?)  and  (3.10).  If 
we  substitute  from  (3.10),  this  gives 


(4,7)  xk£^  =  2  Cijk 


+  2  eij£ 


3g 


VP  7x7  dV 


sjk(5>  If-  d'-<» 


where  we  used  the  Green-Gauss  theorem  and  set  a  surface  term  at  infinity 
to  zero. 

From  (4.7),  one  can  obtain  various  special  cases  involving  surface 
and  line  distributions  of  dislocations.  For  example,  for  a  line  distribu¬ 
tion  of  dislocation  along  a  closed  curve  C  ,  we  obtain 


(4.8) 


'ki 


=  2  Ckijl3j  ^  ^  +  7  (J) 


*i 


2  Jlijuj 


9x.  k 

l 


where  is  the  Burger's  vector  per  unit  length  of  C  and  dj^  is 

the  element  of  the  arc. 

Upon  substituting  (4.8)  into  (3.14),  we  obtain  the  stress  field  due 
to  a  line  distribution  of  dislocations 


(4.9)  tk£/2u  =  Kijbjt 


*  <c,kti  -  VdV 


This  result  is  identical  to  the  Peach  and  Koehler7  formula  with  modifi¬ 
cation  that  here  G  is  the  nonlocal  Green's  function  (4.4)  with  e^O. 
As  we  shall  see,  the  most  interesting  new  feature  of  (4.9)  is  that,  with 
G  given  by  (4.4),  at  a  point  on  the  dislocation  line  C  ,  the  stress  is 
finite  so  that  the  self-stress  and  energies  of  dislocation  loops  can  be 
calculated,  free  of  infinities. 


(ii)  Two-Dimensional  Infinite  Plane 

In  the  case  of  the  plane  strain,  introducing  the  Airy's  stress 
function  ^x^Xj)  by 

*11  *  *,22  ’  *22  “  *,ll  *  t12  *,12 


(4.10) 


12 


we  obtain  an  equation  replacing  (3.16) 
(4.11)  (l  -  e2v2)  V4$  =  2un 


where 


n  =  n33  *  a23,l  "  a13,2 

a23  =  a21 ,2  '  a22,l  ’  ai3  =  all,2  "  ai2,l 


depend  on  x1  and  x2  only. 

Green's  function  in  this  case,  can  be  found  similar  to  decomposition 
(4.2)  with  V2  given  by 


Hence, 


G0  (x-C)-(x-Q 


(4.13)  G(  |  x-£  | )  »  2%  K0(|x-g|/E)  -  —  2  -  -  £n(|x-p/e),  e  t  0 

8ire 


(4.14)  G( |x-C  | ) 


(x-p  •  (x-p 
8tt 


£n(|x-p)  ,  c  =  0 


where  Kq(2)  is  the  modified  Bessel's  function.  Again,  we  take  GQ = 1 
to  render  t^  regular  at  x*  £  . 


( 


l 


V 


Airy's  stress  function  is  obtained  to  be 


(4.15) 


*(x) 


’ 

2U 

. 


S 


tG,l  a23«>  -  G,2  ds 


where  we  used  the  Green-Gauss  theorem  and  set  a  line  integral  to  zero  at 
infinity.  For  a  line  distribution  of  dislocations  in  the  =  0  -  plane, 
we  obtain 


(4.16) 


«(x)  =  -  2p 


[C  1  b2(p  dCj  ♦  G>2  b1(C)  dC2l 


The  stress  field  follows  from  (4.10) 


(4.17) 


tll  *  *  2y 


X22  =  '  2vJ 


*12  s  2y 


(G,122  b2  d^l  +  G,222  bl  dC2}  ’ 
(G,111  b2  d^l  +  G,211  bl  ’ 


(G,112  b2  d?l  +  G,212  bl  d^2^ 


(iii)  Anti-Plane  Strain 

In  the  case  of  anti-plane  strain,  equations  of  equilibrium  are 
satisfied  if 


(4.18) 


*13 


f 


and  we  obtain 


(4.19) 


(1  -  e272)  V2<$>  =  u  a.3 


where 


(4.20) 


a33  =  a31,2  "  a32,l 


Green's  function  for  this  case,  is  obtained  to  be 


(4.21)  G(|x-§|)  =  -  ±  [*n(|x-£|/e)  ♦  K0((x-5|/e)]  ,  e  t  0 


(4.22)  G(|x-5|)  =  -  ±  ^(U-SD 


The  stress  field  is  given  by 


(4.23)  tJ3  =  v 


g>2  b(£)  dcx  d£2 


t23  =  -  y  J  G,1  b(P  dSl  d^2 
S 


For  a  line  distribution  of  dislocations  on  the  plane  x^  =  0  ,  we  have 


(4.24)  t13 


r 

u 

c 


G  -  b(£)  d£  , 
— 


*23 


G  .  b(£)  d£ 
»  1  — 


In  plane  polar  coordinates  (r,0),  the  stress  field  is  given  by 


(4.25) 


Sr  =  u  J  F  If  B<»  dl 

c 


U  b(|)  di 


5.  UNIFORM  DISTRIBUTIONS  OF  SCREWS 

Here  we  calculate  the  stress  field  for  two  different  uniform  line 
distributions  of  screw  dislocations: 


(i)  Screw  Dislocation  Along  a  Straight  Line 

Consider  a  line  distribution  of  screw  dislocations  of  constant 
Burger's  vector  along  a  straight  line  segment  ]  xj<£  ,  x2sx5=0  .  Green’s 

function  is  given  by 


(5.1)  G(|x-§|)  =  -  ^  [In  (p/e)  ♦  KqCd/c)  ] 


where 


(5.2) 


P  =  I  (Xj  -  £)2  +  x2)^ 


We  can  evaluate  t^j  given  by  (4.24)  immediately  since 


=  3G/3x^  =  -  3G/3£ 


and  we  have 


! 
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(5.3)  t23  =  yb  [  G(  k  -  Cl)  -  G(  |2  +  C  i)  ] 


Explicitly, 


(5.4)  t23 


£  {*» 


r  (X1~H)2+X22  -j*S 


L  (x1  +  £)2+X22 


+  K. 


/(x1-£)2+X22/e 


-  K. 


/(x1+^)2+x22/t 


Along  the  line  of  screws,  this  gives 


(s.s)  t23(xvo)  -  -  tn  (-[^r > ♦  ko<!v£!/e> 


Calculation  of  t^3  is  conylicated.  However,  it  is  easy  to  see  that 

(5.6)  t13(x1 ,0)  *  0 

Unlike  the  classical  case,  *23^*1 has  no  singularity  at  the  end 
points  Xj *  ± l  of  the  screw  line.  In  fact,  we  have 

(5.7)  t23(±£,0)  *  ±H_[£n(£/e)  +  KQ(Zi/t)] 
for  l  't»  1  ,  we  have  the  asymptotic  value 


(5.8) 


t23(±Jl,0)  +  ±  ^  [£n(£/e)  ♦  (7re/4i ) 15  exp(-2£/e)J 
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where  the  second  term  can  also  be  neglected  as  compared  to  the  first  one 
for  large  9/  z. 

(ii)  Uniform  Distribution  of  Screw  Disloaations  Along  a  Circle 

Suppose  that  in  the  plane  x^ = 0  ,  there  is  a  uniformly  distributed 
screw  dislocations  along  a  circle  of  radius  R.  In  plane  polar  coordi¬ 
nates,  we  have 

(5.9)  Xj  =  r  cos  0  ,  x^  =  r  sin  0  , 

Cj  !  R  cos  ,  C2  =  R  sin  P  , 

| x  -  C|  =  [  r2  +  R2  -  2r  R  cos  (p-6)  ]'2 


Greens'  function  is  given  by 

(5.1°)  G  (|x  -  §1)  •  -  2W~  t£n(l?  '  +  k0(l?  *  f l/e)] 

To  calculate  the  stress  field  t  and  t  -  we  must  evaluate  two  integrals 

zr  z6  6 

(5*H)  Ij  *  S.n  |  [r^  +  R^  -  2  r  R  cos  ( cp- 0 ) ] ^  |  dp 

6 

27T+S 

(5.12)  12  ~  {  K0  j  J  [r2  +  R"  -  2  r  R  cos  (p-0)  J*5  j  dp 

e 
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To  evaluate  (5.11),  we  write 

C5.133  2n  {  t  [  r2  +  R2  -  2  r  R  cos  (^-e)]1*  j  =  £n (  -Lr-*-R  1 ) 

•  +  iln  (l-aQ  cos2  j)*5 

where 

(5.14)  aQ  =  4  r  R/(r  +  R)2  ,  <fi  -  Q  =  <p 


The  second  term  in  (5.13)  can  be  integrated  by  writing  x  =  cos('4/2)  and 
consulting  Ref.  9,  p.  562,  No.  38.  Consequently 


(5.15) 


2tt  2n(r/e) 
2tt  «n(R/c) 


r  >  R 
r  <  R 


To  evaluate  ,  we  employ  the  integrals  6.684  on  p.  741  of  Ref.  9,  and 
note  that 


where  JQ 


KqU) 

and  Nq 


=  —  [  JQ  (it)  ♦  i  Nq  (it)  ] 

are  zero  order  Bessel  functions. 


The  result  is 


(5.16) 


2tt  Io(R/c)  Ko(r/e) 
2tt  lQ(r/e)  KQ(R/e) 


r  >  R 
r  <  R 


Consequently, 


( 


(5.17)  1  =  J  GUj'§!  dl  =  ' 


£n  (r/e)  +  IQ (R/e)  K^Cr/e),  r>R 
in  (R/e)  +  IQ(r/e)  Ko(R/e),  K* 


The  stress  field  is  given  by 


(5.18) 


bR  31  . 

Sr  =  y  “  *36  =  °  * 


t  „  =  pbR 
26 


~=  { 
3r  \ 


p  b  |  [  f  ”  I0  (R/E)  S(r/e)  ' 

y  b  |  I  (r/e)  K0(R/e)  , 


In  the  special  case  when  R-*0  and  2^Rb  =  bn  ,  we  obtain 


(5.19) 


Se  =  2nr  [  1  '  e  Kl(r/£)  J 


which  is  identical  to  our  previous  result  for  a  single  screw  dislocation 
having  Burger’s  fector  bn  . 


-•ms; 


6.  STRESS  DISTRIBUTIONS 


Here,  I  present  some  numerical  results  on  the  stress  distributions 
for  the  cases  discussed  in  Section  5  and  establishes  a  fracture  criteria 
based  on  the  maximum  shear  stress* 

(i)  Single  Screw:  The  shear  stTess  given  by  (5.19)  may  be  expressed 
in  non-dimensional  form: 


(6.1) 


-1 , 


T0(p)  *  (2re/yb)  t,Q  =  P  X[1  -OKj (p) ] 


where 

(6.2) 


p  =  r/e 


The  stress  field  given  by  (b.l)  is  displayed  graphically  in  Fig.  1.  It  has 
no  singularity  at  p  =  0.  In  fact,  TQ(p)  vanishes  at  p-0  in  contradic- 
tion  to  the  classical  elasticity  solution  which  gives  infinite  stress  at  c=  0 
The  maximum  stress  occurs  at  p  =  1.1  and  is  given  by 


(6.3) 


:9  max 


■  »• 3993  nr 


If  we  write  h= e/0.3993,  this  agrees  with  Frenkel's  estimate  of  the  theore¬ 
tical  strength  of  single  crystals,  based  on  atomic  considerations  (cf., 

Kelly  [10],  p.  12).  In  fact,  if  we  use  e  *  ega  =  0.39  a  ,  which  is  obtained 
on  the  basis  of  matching  of  the  dispersion  curve  predicted  by  non-local  elas 
ticity  and  the  Bom-Karman  lattice  model^,  we  find  for  the  single  aluminum 


crystal 
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(6.4)  tC/u  =  0.12  (Al:  £111]  <  110  >  } 

This  is  very  close  to  the  theoretical  strength  t  /p  =  0.11  based  on 
atomic  models. 


(ii)  Screw  Dislocations  Along  a  Straight  Line  Segment:  Even  single 
crystals  contain  many  dislocations.  For  a  uniform  distribution  of  screws 
along  a  line  segment  | Xj |  < £  ,  x,  =  x,  = 0  ,  the  shear  stress  given  by 

(5.5)  may  be  written  in  non-dimensional  form 


(o.5) 


f2  =  t25/td 


in 


I  x+1  1 


|X-1  i 


Kn  (y  |  x+1  j )  -  K  (y|x-r) 


where 

(6.6)  t  =  ub/2-  =  ub0N7-ti  ,  x  =  Xj/i  ,  y  -  Uc 


Here  bQ  is  the  atomic  Burger's  vector  and  N  is  the  total  number  of 
dislocations  over  a  distance  £  . 

The  distribution  of  the  shear  stress  (6.5)  as  a  function  of  x  is 
shown  in  Fig.  2  for  various  values  of  y  .  Behavior  of  T-,  is  governed 
basically  by  the  first  term  in  (6.5)  except  near  x  =  l  .  At  x  =  l  ,  we 
have 


(6.7) 


t,i(!) 


U  b  N 
o 

2~e 


In  y 


The  value  of  T^(l) 


is  very  close  to  the  maximum  stress  for 


>  >  3 


i 


(cf.  Table  1).  For  a  single  atomic  dislocation  according  to  (6,3),  we 
have  the  theoretical  strength 


(6.8) 


c  H  0 

t  =  0.3993  ~ 

v  2te 


Combining  (6.7)  and  (6.8),  we  obtain 


td/tC 

y  y 


N  In  y 
0.3993  y 


where  we  set  t-,-(l)  =  t  =  the  yield  stress  for  the  distributed  dis¬ 
locations.  This  gives  the  shear  stress  reduction  due  to  the  presence 

of  2N  dislocations  distributed  uniformly  along  a  straight  line  segment 

d  c 

of  length  2i  .  Since  ty  <_  t^  ,  the  maximum  number  of  dislocations  is 
given  by 


(6.10) 


N  =0. 3995  - - 

max  cn  > 


For  v  =4.02  *  10  ,  this  gives  N  =  1514  ,  which  mav  be  conservative 

°  max 

since  the  distribution  is  not  generally  uniform  but  in  an  inverse  pile-up 
configuration^ . 


(iii)  Uniform  Distribution  of  Screws  Along  a  Circle:  In  this  case, 
the  stress  fields  given  by  (5.18)  may  be  expressed  in  non-dimensional  form 


10.  ID 


=  t.p/ub  = 


-  -  <  I0(<)  ^ (<w)  , 


T 


k  K0(<)  Ij(<C) 


o  <  1 
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where 

(6.12)  P  =  r/R  ,  k  =  R/c 


T  as  a  function  of  p  ,  for  various  values  of  <  ,  is  displayed  in 

Fig.  3.  For  value  of  <  50,  the  maxima  of  T  occurs  near  p  =  1  . 

The  locations  and  values  of  T  are  given  in  Table  2. 

max 

If  bQ  is  the  atomic  Burger's  vector,  then 


(6.15)  2rr  R  b  =  N  bQ 

where  N  is  the  number  of  dislocations  on  the  dislocation  circle  with 
radius  R.  Using  (6.8),  (6.12)  and  (6.13),  we  obtain 


(6.14) 


tz6/ty 


N 

0.3993  < 


T 


According  to  Fig.  3,  0.324  <_  T^ax  £  1  .  Consequently, 

(6.15)  N / 1 • 2  <  <  ty/ty  <  N/0.4  < 


For  perfect  crystals  with  e  =0.39  a  ,  this  gives  approximately 

(6.16)  0.5  Na/R  <  t^/tC  <  Na/R 

y  y  - 


indicating  reduction  of  the  yield  stress  with  the  presence  of  large 
number  dislocations  uniformly  distributed  along  a  circle  of  radius  R  . 


Table  1  :  Maximum  Shear  Stress  and  its  Location 


(Line  Segment) 


1 

1.5 

2 

3 

5 

10 

1.446 

1.197 

1.103 

1.039 

1.000 

1.000 

0.7478 

1.0501 

1.3008 

1.6851 

2.3026 

2.9957 

Table  2  :  Maximum  Shear  Stress  and  its  Location 


(Circle) 


<  - 

1 

2 

3 

5 

10 

50 

P  * 

1.8 

1.5 

1.4 

1.3 

1.2 

1.1 

T 

0.3243 

0.4836 

0.5688 

0.6650 

0.7688 

0.9058 

max 
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Universities  (Con't.) 

Professor  Valctr  D.  Pilkey 
University  of  Virginia 
Research  Laboratories  for  the 
Engineering  Sciences  and 
Applied  Sciences 
Charlottesville,  Virginia  22901 

Professor  E.  D.  Wil Inert 
Clarkson  Collage  of  Technology 
Department  of  Mechanical  Engineering 
Potsdam,  Raw  York  13676 

Dr.  Walter  E.  Haisler 
Texas  46 M  University 
Aerospace  Engineering  Department 
College  Station,  Texas  77843 

Dr.  Hussein  A.  Eaael 
University  of  Arixona 
Department  of  Aerospace  and 
Mechanical  Engineering 
Tucson,  Arizona  85721 

Dr.  S.  J.  Fenves 
Ca me gi e-Mel  Ion  University 
Department  of  Civil  Engineering 
Schenley  Park 

Pittsburgh,  Pennsylvania  15213 

Dr.  Ronald  L.  Huston 
Department  of  Engineering  Analysis 
University  of  Cincinnati 
Cincinnati,  Ohio  45221 

Professor  C.  C.  M.  Sih 
Lehigh  University 
Institute  of  Fracture  and 
Solid  Mechanics 
Bethlehem,  Pennsylvania  18015 

Professor  Albert  S.  Kobayashi 
University  of  Washington 
Department  of  Mechanical  Engineering 
Seattle,  Washington  98105 

Professor  Daniel  Frederick 
Virginia  Polytechnic  Institute  and 
State  University 

Department  of  Engineering  Mechanics 
Blacksburg,  Virginia  24061 


Universities  (Con't) 

Professor  A.  C.  Eringen 
Princeton  University 
Department  of  Aerospace  and 
Mechanical  Sciences 
Princeton,  Rev  Jersey  08540 

Professor  E.  B.  Lee 

Stanford  University 

Division  of  Engineering  Mechanics 

Stanford,  California  94305 

Professor  Albert  I.  Ring 
Wayne  State  University 
Biomechanics  Research  Center 
Detroit,  Michigan  48202 

Dr.  V.  R.  Hodgson 
Wayne  State  University 
School  of  Medicine 
Detroit,  Michigan  48202 

Dean  B.  A.  Boley 
Northwestern  University 
Department  of  Civil  Engineering 
Evanston,  Illinois  60201 

Professor  P.  C.  Hodge,  Jr. 
University  of  Minnesota 
Department  of  Aerospace  Engineering 
and  Mechanics 

Minneapolis,  Minnesota  55455 

Dr.  D.  C.  Drucker 
University  of  Illinois 
Dean  of  Engineering 
Urbane,  Illinois  61801 

Professor  R.  M.  Nevmark 
University  of  Illinois 
Department  of  Civil  Engineering 
Urbane,  Illinois  61803 

Professor  E.  Raissner 
University  of  California,  San  Diego 
Department  of  Applied  Mechanics 
La  Jolla,  California  92037 

Professor  Willima  A.  Rash 
University  of  Massachusetts 
Department  of  Mechanics  and 
Aerospace  Engineering 
Bmhkrst,  Massachusetts  01002 


Pniveraitiea  (Coo'c) 

Profeaaor  C.  Hermann 
Stanford  Oniveraity 
Department  of  Applied  Mechanic a 
Stanford,  California  94305 

Profeaaor  J.  D.  Achenbach 
Northweat  Oniveraity 
Department  of  Civil  Engineering 
Evanaton,  lllinoia  60201 

Profeaaor  S.  B.  Dong 
Oniveraity  of  California 
Department  of  Mechanica 
Loa  Angelea,  California  90024 

Profeaaor  Burt  Paul 
Oniveraity  of  Pennaylvania 
Tovne  School  of  Civil  and 
Mechanical  Engineering 
Philadelphia,  Pennaylvania  19104 

Profeaaor  H.  W.  Liu 
Syracuee  Oniveraity 
Department  of  Chemical  Engineering 
and  Metallurgy 
Syracuae,  New  York  13210 

Profeaaor  S.  Bodner 
Technion  BAD  Foundation 
Haifa,  larael 

Profeaaor  Werner  Coldamith 
Oniveraity  of  California 
Department  of  Mechanical  Engineering 
Berkeley,  California  94720 

Profeaaor  R.  S.  Rivlin 
Lehigh  Oniveraity 
Center  for  the  Application 
of  Mathemat ica 

Bethlehem,  Pennaylvania  18015 

Profeaaor  F.  A.  Cozzarelli 
State  Oniveraity  of  Mew  York  at 
Buffalo 

Diviaion  of  Interdiaciplinary  Studiea 
Harr  Parker  Engineering  Building 
Chemiatry  Road 
Buffalo,  New  York  14214 
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Oniveraitiea  (Con't) 

Profeaaor  Joaeph  L.  Roae 
Drexe 1  Oniveraity 

Department  of  Mechanical  Engineering 
and  Mechanica 

Philadelphia,  Pennaylvania  19104 

Profeaaor  B.  K.  Donaldaon 
Oniveraity  of  Maryland 
Aeroapace  Engineering  Department 
College  Park,  Maryland  20742 

Profeaaor  Joaeph  A.  Clark 
Catholic  Oniveraity  of  America 
Department  of  Mechanical  Engineering 
Waahington,  D.C.  20064 

Dr.  Samuel  B.  Bat dor f 
Oniveraity  of  California 
School  of  Engineering 
and  Applied  Science 
Loa  Angelea,  California  90024 

Profeaaor  laaac  Fried 
Boaton  Oniveraity 
Department  of  Mathemat ica 
Boaton,  Maaaachuaetta  02215 

Profeaaor  E.  Krempl 
Rensaelaer  Polytechnic  Inatitute 
Diviaion  of  Engineering 
Engineering  Mechanica 
Troy,  Hew  York  12181 

Dr.  Jack  R.  Vinaon 
Oniveraity  of  Delaware 
Department  of  Mechanical  and  Aeroapace 
Engineering  and  the  Center  for 
Compoaite  Materiala 
Newark,  Delaware  19711 

Dr.  J.  Duffy 
Brown  Oniveraity 
Diviaion  of  Engineering 
Providence,  Rhode  Ialand  02912 

Dr.  J.  L.  I wed low 
Carnegie-Mellon  Oniveraity 
Department  of  Mechanical  Engineering 
Pittaburgh,  Pennaylvania  15213 
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Universities  (Con't) 

Dr.  V.  K.  Varadan 

Ohio  State  University  Reaearch  Foundation 
Department  of  Engineering  Mechanics 
Columbus,  Ohio  43210 

Dr.  Z.  Haahin 

Dniveraity  of  Pennaylvania 
Department  of  Metallurgy  and 
Materiala  Science 
College  of  Engineering  and 
Applied  Science 

Philadelphia,  Pennsylvania  19104 

Dr.  Jackson  C.  S.  Tang 
University  of  Maryland 
Department  of  Mechanical  Engineering 
College  Park,  Maryland  20742 

Professor  T.  T.  Chang 
University  of  Akron 
Department  of  Civil  Engineering 
Akron,  Ohio  44325 

Professor  Charles  V.  Bert 
University  of  Oklahoma 
School  of  Aerospace,  Mechanical, 
and  Nuclear  Engineering 
Norman,  Oklahoma  73019 

Professor  Satya  N.  Atluri 
Georgia  Institute  of  Technology 
School  of  Engineering  and 
Mechanics 

Atlanta,  Georgia  30332 

Professor  Graham  F.  Carey 
University  of  Texas  at  Austin 
Department  of  Aerospace  Engineering 
and  Engineering  Mechanics 
Austin,  Texas  78712 

Dr.  S.  S.  Vang 
University  of  Illinois 
Deoartment  of  Theoretical  and 
Applied  Mechanics 
Urbana,  Illinois  61801 

Professor  J.  F.  Abel 
Cornell  University 
Department  of  Theoretical 
and  Applied  Mechanics 
Ithaca,  Nev  Tork  14853 


Universit ies  (Con't) 

Professor  V.  B.  Neubert 
Pennsylvania  State  University 
Department  of  Engineering  Science 
and  Mechanics 

University  Park,  Pennsylvania  16802 

Professor  A.  V.  Leiasa 

Ohio  State  University 

Department  of  Engineering  Mechanics 

Columbus,  Ohio  43212 

Professor  C.  A.  Brebbia 
University  of  California,  Irvine 
Department  of  Civil  Engineering 
School  of  Engineering 
Irvine,  California  92717 

Dr.  Ceorge  T.  Hahn 
Vanderbilt  University 
Mechanical  Engineering  and 
Materials  Science 
Nashville,  Tennessee  37235 

Dean  Richard  B.  Gallagher 
University  of  Arizona 
College  of  Engineering 
Tucson,  Arizona  85721 

Professor  E.  F.  Rybicki 
The  University  of  Tulsa 
Department  of  Mechanical  Engineering 
Tulsa,  Oklahoma  74104 

Dr.  R.  Raftka 

Illinois  Institute  of  Technology 
Department  of  Mechanics  and  Mechanical 
and  Aerospace  Engineering 
Chicago,  Illinois  60616 

Professor  J.  G.  de  Oliveira 
Massachusetts  Institute  of  Technology 
Department  of  Ocean  Engineering 
77  Massachusetts  Avenue 
Cambridge,  Massachusetts  02139 

Dr.  Bernard  V.  Shaffer 
Polytechnic  Institute  of  New  York 
Route  110 

Farmingdale,  New  Tork  11735 
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Industry  and  Research  Institutes 

Dr.  Norman  Hobbs 
Reman  AviDyne 
Division  of  Raman 
Sciences  Corporation 
Burlington,  Massachusetts  01803 

Argonne  National  Laboratory 
Library  Services  Department 
9700  South  Caas  Avenue 
Argonne,  Illinois  60440 

Dr.  M.  C.  Junger 
Cambridge  Acoustical  Associates 
54  Rindge  Avenue  Extension 
Cambridge,  Massachusetts  02140 

Mr.  J.  R.  Torrance 
General  Dynamics  Corporation 
Electric  Boat  Division 
Groton,  Connecticut  06340 

Dr.  J.  E.  Greenspon 

J.  G.  Engineering  Research  Associates 
3831  Menlo  Drive 
Baltimore,  Maryland  21215 

Newport  News  Shipbuilding  and 
Dry  Dock  Company 
Library 

Newport  News,  Virginia  23607 

Dr.  W.  F.  Bozich 

McDonnell  Douglas  Corporation 

5301  Bolsa  Avenue 

Huntington  Beach,  California  92647 

Dr.  B.  N.  Abramson 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio,  Texas  78284 

Dr.  I.  C.  DeHart 
Southwest  Research  Institute 
8500  Culebra  load 
San  Antonio,  Texas  78284 

Dr.  M.  L.  Baron 
Weidlinger  Associates 
110  East  59th  Street 
New  Tork,  New  Tork  10022 


Industry  and  Research  Institutes  (Con' 
Dr.  T.  L.  Geers 

Lockheed  Missiles  and  Space  Company 

3251  Hanover  Street 

Palo  Alto,  California  94304 

Mr.  Villiam  Caywood 
Applied  Physics  Laboratory 
Johns  Hopkins  Road 
Laurel,  Maryland  20610 

Dr.  Robert  E.  Dunham 
Pacifica  Technology 
P.0.  Box  148 

Del  Mar,  California  92014 

Dr.  M.  P.  Kanninen 
Battelle  Columbus  Laboratories 
505  King  Avenue 
Columbus,  Ohio  43201 

Dr.  A.  A.  Hochrein 
Daedalean  Associates,  Inc. 

Springlake  Research  Road 
15110  Frederick  Road 
Woodbine,  Maryland  21797 

Dr.  James  W.  Jones 
Swanson  Service  Corporation 
P.0.  Box  5415 

Huntington  Beach,  California  92646 

Dr.  Robert  E.  Nickell 
Applied  Science  and  Technology 
3344  North  TorTey  Pines  Court 
Suite  220 

La  Jolla,  California  92037 

Dr.  Kevin  Thomas 
Westinghous%  Electric  Corp. 

Advanced  Reactors  Division 
P.  0.  Box  158 

Madison,  Pennsylvania  15663 

Dr.  H.  D.  Bibbitt 
Hibbitt  6  Karlsson,  Inc. 

132  George  M.  Cohan  Boulevard 
Providence,  Rhode  Island  02903 

Dr.  R.  D.  Mindlin 
89  Deer  Hill  Drive 
Ridgefield,  Connecticut'  06877 
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Or.  Norman  Hobbs 
Inua  AviDyne 
Division  of  Raman 
Sciences  Corporation 
Burlington,  Massachusetts  01803 

Argonne  National  Laboratory 
Library  Services  Department 
9700  South  Cass  Avenue 
Argonne,  Illinois  60440 

Dr.  M.  C.  Junger 
Cambridge  Acoustical  Associates 
54  Rindge  Avenue  Extension 
Cambridge,  Massachusetts  02140 

Mr.  J.  H.  Torrance 
Ceneral  Dynamics  Corporation 
Electric  Boat  Division 
Croton,  Connecticut  06340 

Dr.  J.  E.  Greenspon 
J.  C.  Engineering  Research  Associates 
3831  Menlo  Drive 
Baltimore,  Maryland  21215 

Newport  News  Shipbuilding  and 
Dry  Dock  Company 
Library 

Newport  News,  Virginia  23607 

Dr.  W.  F.  Bozich 

McDonnell  Douglas  Corporation 

5301  Bolsa  Avenue 

Huntington  Beach,  California  92647 

Dr.  H.  N.  Abramson 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio,  Texas  78284 

Dr.  R.  C.  DeHart 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio,  Texas  78284 

Dr.  M.  L.  Baron 
Weidlinger  Associates 
'10  East  59th  Street 
New  Tork,  New  York  10022 
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Industry  and  Research  Institutes  (Con't) 


Dr.  T.  L.  Geers 

Lockheed  Missiles  and  Space  Company 

3251  Hanover  Street 

Palo  Alto,  California  94304 

Mr.  William  Caywood 
Applied  Physics  Laboratory 
Johns  Hopkins  Road 
Laurel,  Maryland  20810 

Dr.  Robert  E.  Dunham 
Pacifica  Technology 
P.0.  Box  148 

Del  Mar,  California  92014 

Dr.  M.  F.  Eanninen 
Battelle  Columbus  Laboratories 
505  Ring  Avenue 
Columbus,  Ohio  43201 

Dr.  A.  A.  Hochrein 
Daedalean  Associates,  Inc. 
Springlake  Research  Road 
15110  Frederick  Road 
Woodbine,  Maryland  21797 

Dr.  James  W.  Jones 
Swanson  Service  Corporation 
P.0.  Box  5415 

Huntington  Beach,  California  92646 

Dr.  Robert  E.  Hickell 
Applied  Science  and  Technology 
3344  North  Torrey  Pines  Court 
Suite  220 

La  Jolla,  California  92037 

Dr.  Kevin  Thornes 
Westinghouse  Electric  Corp. 

Advanced  Reactors  Division 
P.  0.  Box  158 

Madison,  Pennsylvania  15663 

Dr.  B.  D.  libbitt 
Hibbitt  6  Ear Is  son,  Inc. 

132  George  M.  Cohan  Boulevard 
Providence,  Rhode  Island  02903 

Dr.  R.  D.  Hindi in 
89  Deer  Hill  Drive 
Ridgefield,  Connecticut  06877 
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Induatrr  «ad  Reeearch  Inatitute ,  (Con't) 

Dr.  Richard  E.  Daae 
Mega  Engineering 
11961  Tech  Road 

Silver  Spring,  Maryland  20904 

Mr.  C.  M.  Stanley 
Lockheed  Palo  Alto  Reaearch 
Laboratory 
3251  Hanover  Street 
Palo  Alto,  California  94304 

Mr .  R.  t, .  Cloud 

®°kert  L.  Cloud  Aaeociatea,  Inc. 

2972  Adeline  Street 
Berkeley,  California  94703 


